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Abstract 

It has long been known how to construct radiation boundary conditions 
for the time dependent wave equation. Although arguments suggesting 
that they are accurate have been given, it is only recently that rigorous 
error bounds have been proved. Previous estimates show that the error 
caused by these methods behaves like eC^e"^*! for any 7 > 0. We improve 
these results and show that the error behaves like Cet^ . 

1 Introduction 

Numerical solution of time dependent wave equations is an important problem 
in physics, engineering and mathematics. To solve the wave equation on M'^+^, 
one must truncate the domain to a finite region due to the limited memory 
of most computers. Of course, on a finite region, boundary conditions must be 
specified in such a way as to minimize spurious reflections. Boundary conditions 
of this form were first described in [Ij [2l S [H [5] , although rigorous error bounds 
would wait until more recently [51 [7] . 

In [3] , a family of absorbing boundary conditions based on rational function 
approximation to the Dirichlet-to-Neumann operator in the frequency domain 
are reviewed. Boundary conditions for the half-space (a boundary at a; = 0), as 
well as cylindrical and spherical coordinates are also constructed. Error bounds 
are proved for this family by inverting the Fourier-Laplace transform for both 
the true solution and his approximation and bounding the difference. Due to 
poles on the imaginary line in s (s being the variable dual to t), the difference is 
bounded on a countour separated from the singular points, namely a line in the 
right half plane 7 -I- iM. This shows that the error is bounded by C^e^*, with 

left implicit. 

A careful examination of the poles of the rational function reveal that they 
approximate the branch cut of the true solution in the sense of hyperfunctions 
[S]. Instead of using the machinery of hyperfunctions, we take an elementary 
approach. The true solution can be represented as a certain integral over a com- 
pact region. The approximate solution, after we collect the residues associated 
to the poles on the imaginary line, turns out to be a quadrature for this integral. 
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By computing the difference between tfie quadrature and the true integral, we 
can compute an optimal error bound. 

Let us now state our results precisely. Let u{x, y, t) solve: 

dtu{x, y, t) = dlu{x, y, t) + l^yu{x, y, t) (1) 

where a; € M and y £ R^^^ (x is the normal direction, y the tangential direc- 
tions). 

We wish to solve ([T]) on M^+^. The boundary will be taken to be the surface 
1 = 0, and thus the approximation region will be the region {(x, y, z,t) : x > 0}. 
We let ui,{x,y,t) be the approximation, solving ([T]) on the half-space. The 
boundary conditions imposed are Hagstrom's: 

n (^°^ (^) ' = ° 

The main theorem is the following: 

Theorem 1. Let u{x,y,t) solve ^ on and Ub{x,y,t) solve ([T]) with 

boundary conditions Q . Then we have the following error bound: 

\u{x,y,t) - Ub{x,y,t)\ 

< ^ ^ (2^) (2"*' + ^rit + nKmax 8n -h 3) 

K 



^0\^-^{K^a.+i')j (3) 

2 Proof 

2.1 The Exact Boundary 

We begin by reviewing the exact boundary conditions described in [B] . Applying 
the Laplace transform of ([T]) with respect to time (letting s be dual to t) and 
the Fourier transform with respect to y (with k dual to y) yields: 

s^Ji = dlii - k^t (4) 

The solution to (g]) is: 

Ti{x) = A{s, k)eV-'+|k|^ + B{s, k)e-V''' + |kP^ (5) 

The solutions with nonzero A(s,k) are nonphysical, since they correspond to a 
wave coming from infinity to the object. Thus our boundary conditions must 
imply A{s,\i) — 0. Such a boundary condition is (in the frequency domain): 



d^f[{x, k, s) + v/.s2 + |k|2t(a;, k, s) = (6) 
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Of course, the operator + |kp is non-local in time and space, so we will 
approximate it. 

To reduce the dependence to a single variable, we make the substitution 
z = s/ |k|, yielding: 

dxu{x, k, s) + |k|-\/l + z^u{x, k, s) = 
This boundary condition can be rewritten as: 

dx%{x, k, s) + |k| (z + ^=^) Hx, k, s) = (7) 

Let h{z) = |k| /{z + y/l + z^). We will invert the Laplace transform, and 
shift the contour to surroung the singularities of h{z). Ths following lemma 
summarizes the necessary analyticity properties of h(z): 

Lemma 1. The function h{z) is analytic on<C \ [~i,i]. In addition, the differ- 
ence across the branch cut is given by: 

lim {h{z + e) - h{z - e)) = 2|k| ^1 + z^ (8) 

Proof. The function h{z) is well defined and analytic for iRz > 0. It is strictly 
imaginary on {z : Jftz = Oand|z| > 1}. By the Schwartz reflection principle, it 
can be analytically continued to the left half plane, with a discontinuity along 
the line [—i, i]. 

An explicit calculation shows ([5]). □ 

We now reconstruct u{x,y,t). This is done by inverting the Laplace trans- 
form: 

uix,y,t) = ,fl,^.,^ [ e^* / e'y^ll{x,k,s)dkds (9a) 
?l(x,k, t) = / e''*fi{x,'k,s)ds (9b) 

J a+iR 

And so, the integral we must approximate is 

/ 2|kiym2/(z)e^*dz. 

J —i 



2.2 The Approximation 

We review the approximation itself, and how ([2]) was derived. Our description 
follows [B] quite closely. We approximate |k| Vl + 2^ by: 



HvTT^Hk|(, + ^:p^^j«=|k|| = + 5^^| (10) 

' ■2z ■ 
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where the right hand side is the n'th iteration of the continued fraction. 
A straightforward computation shows that in the frequency domain, 

/ ^ 

z H — I u{x, k, s) = 



9a;u(a;,k,s) + |k| , 

2z i 

V 

corresponds to the boundary condition We simphfy this 



Lemma 2. Let 9j = jii / {n + 1). Then we have the following formula: 

1 " -2/1 

1 ^ sm^, ^^^^ 
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Proof. Let Un{x) be the n*'' Chebyshev Polynomial of the 2"'^ kind and P„(z) be 
the successive numerators of the sequence of finite continued fractions {Po{z) = 
l,Pi(z) = 2z,...). 

If we take Un{iz) then for n odd we get iPn{z) and for n even we get Pn{z), 
and the sequence of finite continued fractions is (P„_i(z))/(P„(z)) for n > 1. 

We consider the case case where n is even; in this case, the finite continued 
fraction is {iPn-i{z)){Pn{z))- Thus, ratios of Chebychev polynomials of the 2"'^ 
kind only differ from the finite continued fraction by multiplication by i, and 
will therefore have the same zeros. The continued fraction will have poles where 
P„(z) is zero, for n even. That is, when C/„(iz) = 0. 

We will take z = i cos 6. Thus, we are looking for zeroes of [/„(— cos 6) where 
C/„_i(-cos0) ^0. 



[/„ (- cos 9) = [/„ (cos 6)^ Unix) ^ smin + l)0 ^ g ^ q^^^2tt, . . . 

smO 

So sin(rt + 1)9 = are our solutions, and that is 6* = — 9j as claimed. 

Thus, U {i cos 9j) =0 and hence z — i cos 9j are the only poles of the contin- 
ued fraction approximation. A direct computation shows that the residues at 
the pole icosOj is [siri^ 9 j) / [n + 1). □ 



As the continued fraction is a close approximation to vT+z^, we can use it 
to approximate an integral involving \/\ + z"^ by substituting the approximation, 
which is a rational function. And so, in evaluating the integral around the 
branch cut, a finite sum which approximates this integral is given by the sum 
of the residues at the poles of the rational function above. 

2.3 The Error Bound 

First, we make the definition |/c|(7(s/|/c|) = u{s, k). 
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Proposition 1. The following error hound holds. 

/i n 
y/lT^e^*g{z)dz - 27Ti\k\ ^ g{zj)e''^*aj 
-i 



< 



3 {n + iy 



:M{x){2nf + %nt + nK^ax + ^n + i) (12) 



Here, Zj — i cos 6j are the positions of the poles, aj the residue at 6j , K^ax 
the maximal frequency under consideration, M{x) is a pointwise upper bound 
on u, dgU and d^u and n is the order of the continued fraction approximation. 

We will need the following lemma 

Lemma 3. 

i-Ax 

f{x)dx - f{0)Ax 



□ 



Proof. This follows immediately from the Intermediate Value Theorem. 

We will prove Lemma [T] above in several intermediate steps. 
Proposition 2. 



Vl + z'^e''*g{z)dz - 27rifc ^ .g(zj)e'^^*a 



<-|k|A03 3 



E 









max g | 




+ 3 


sG[i|k|,i k cos Ae] ' 






max 







max g 

se[-i\k\ COS Ae,-i|k|] 



se[i|k| cos(ej-Ae),i|k| cos(ej+Ae)] 

2 



r,|k| 



k|" 
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+ itg' 



Ikl'' 



2tg' 



2tg' 



s 
s 



5g' 



r,|k| 



r,|k| 



5^.9(TjT,|k| 



s 



Ikl'' 



2g 



Ikl'' 



(13) 



In this equation Zj are the positions of the poles, aj are the residues at the poles. 



A9 = 2{7^+i) -■ '^'^'^ '^(^) ~ g{icos9)e 



it cos 9 f 
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Proof. We first change variables to z = icosO. Using the fact that 



2z- 



2z+ 



we can approximate Vl + by taking a truncated continued fraction. This 
yields: 



|fc| 



„7r " 

2 / 5(i cos 6l)e** " sin^ 6id6' - 27riVg(i cos 6*^)6** 



cose,sin2 6'j 



n + 1 



(14) 



We define f{z) = g{z)e^^ and — 2{n+i) ^^'^ expand the integral around 
each pole to obtain 



/•AO /'TT 

2i / /(icos6i)sin^6id6' + 2i / f{i cos 9) sin 

^0 Jti-A0 



Tl = |fc| 



f{i cos 9) sin^ 6'd6' — ng{i cos9j) 



sm 



n + 1 



(15) 



To simplify further, we substitute a{9) — f{i cos 9) sin^ 9 cancel terms, and use 
the fact that Jg'^^g a'{9j){9 ~ 9j)d9 ^ to get 



pAO pir " / j-Bj+Ae 

2i a{9)d9 + 2i a{9)d9 + 2iy2\ 

Jo Jtt-AB j^-^ \ Js 



(ai9)-ai9,) 



-AS 



a'{9,){9-9,)) \d9 



(16) 



By the triangle inequality, and the mean value theorem, we have 



(lIll)<2|fc|M \a{9)\d9 



-Ad 



Y{9)\d9 



n 



max ct"{0 
^e[ej~Ae,9j+Ae] 



A9^ (17) 



To deal with the ends of the integral, we substitute a and / back into the inte- 
grals near the endpoints. We then use the fact that J^^ f < max^-gj^ f{x){b~o) 
and |sin6'| < \9\ and |sin(7r — 9)\ < — 9\ to obtain: 



(HZD < 2\k\\ 



max 17(1 cos( 

[O,A0]^ 



A9-' 



max g [I cost 

[7r-Ae,7rl 



A9-' 



-y 



max ct"(0 
i£[ej-Ae,ej+Ae] 



A9^ (18) 
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Upon substitution back and simplification, this becomes 



(Ull) < -|A:|A0M 3 



max g{z) 

cos Af?] 



max g{z) 

-■I COS Ay, — ij 



E 



max e^M5"W^' + 2g"(z)z^+g"(z) + 2VW^' 

[icos(ej-A9),jcos(9j+A6()] \ 



+4tg'{z)z^ + 2tg'{z)+rg{z)z'^ + 2rg{z)z^+rg{z) + 5z^g'{z) + 5g'{z)z + 5z^tg{z) 



+ 5tg{z)z + Az'^g(z) + 2g{z) 



(19) 



We also know that g{z) — g (s/|k|, /c). All the derivatives in (PO]) arc in s/|k| 
And so we get: 



HI) < -|k|A03 3 



max g I Ikl 

!e[i|k|,i|k| cos Ae] \ |k| 



max 5 1 777, Ikl 

se[-i|k| COS Ae,-i|k|] \ |k| 



E 



max e 

se[i|k| cos(ej-A6l),i|k| cos(6lj+A6l)] 



k| 

■2t'5(^>|k| 



4^3' T^,|k| +2t.9'^,|k| +t^, -^,|k| 



(20) 



□ 



Now we can prove the final bound, and complete the proof of the main 
theorem. Once Proposition [T] is proven, this implies the main result by (j9ap and 
(I9b|. 

Proof of Proposition [H 

dsu{s, k) = dskg{s/k, k) = kDi\g[s/k, k) = Dig{^^,k). 

dlJi[s,k) = dl{kg{s/k,k)) = kdlg{s/Kk) - l/kDlg{^^,k) 

So Dig = dsM and D\g = |k|9s?I 
Thus, (fT51) above can be simplified to: 
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3 



max u(s, |k|) 

se[j:|k|,i|k| cos Ae] 



max m(s, |k|) 

se[-i|k| COS A6l,-i|k|] 



E 



max 

se[i|k| cos(0j-Ae),i|k| cos(ej+Ae)] 
2 



ewM|k|92^(.,|k|) - 



2\k\diuis, |k|) U- + Ikia^s, |k|) + 2tdsu{s, |k|) — 



+ 4t5,«(s,|k|) — +2ia,u(s,|k|)+t^Il(s,|k|) — +2t'uis,\k\) [ — 



+ t'u{s, |k|) + 5 — a,u(s, |k|) + 59,u(s, |k| 



m{s,\k\) 



5M(s, |k|)— + 4 — |k|) + 2f[{s, |k|) 



(21) 



If we find the maximum for eacli term independently, we will obtain an upper 
bound for this. Noting that cos0 is monotonic decreasing on [0,7r], we obtain: 



IB) < -A0^ ( 3 



max u (s, |k|) 

se[j|k|,i|k| cos AS] 



+ 3 



max u (s, |k|) 

se[-i|k| COS A6l,-i|k|] 



E 



max e 

se[i|k| cos(ej-A0),i|k| cos(ej+A0)] 



icos(93-A6l)t 



|k|92^(s, |k|)cos4(6ij - Ae*) 



2|k|52^(s, |k|) cos^{9j - Ad) + \k\d^,iiis, |k|) + 2tdsfi{s, \k\)cos'^ 



Ad) 



- 4td,uis, |k|) cos''(e'j - AO) + 2td,u{s, |k|) + t^u{s, |k|) cos(6ij - AO) 

- 2t^Ji{s, |k|)cos2(6'j - A6) + t^Ji{s, |k|) - bicos^{ej - Ae)dsTi{s, |k|) 

+ 5dsil{s, \k\)icos{ej-Ae)-5i cos^{ej-A0)tTiis,\k\) + btTi:{s,\k\)icos{9j~Ae) 



4 cos 



A6')^(s,|k|) + 2^(s, |k|) 



(22) 



We also know that u{s, x, k) is bounded, so let U{x) be an upper bound, let 
U'{x) be an upper bound of dgU and U"{x) and upper bound of d^u. Then let 
M{x) be the maximum of these functions. Now, things simplify further to: 
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2 

(m < -A03^/(a;) 6+^ 
V 1=1 

+ 2t cos^{9j - A9) - U cos^ {e.j -M) + 2t + t^ cos(6ij -M)- 2f cos^ {Oj -Ae) + f 
- bicoa^iOj - AO) + 5icos(6ij - AO) - bicos^iOj - Ae)t + 5ti coaie^ - AO) 



4cos^{9j - Ae) + 2 



(23) 



Applying the triangle inequality, we obtain the following as a bound. 



(ESI) < -Ae^M{x) ( 6 + ^ (|k| cos^(0j - A0) ~ 2|k| co-&^{Bj - A0) + |k| 

+ 2t cos^ {9 J - AO) - 4< cos^ (6ij -Ae) + 2t + t^\ cos{9j -Ae)\- 2f cos^ (Oj - AO) 
+ f -5\ cos^iOj -Ae)\ + 5| cos(6ij - A6i)| - 5| 008^(61^ - Ae)\t 

+ 5t\ cos{0j - Ae) \ - 4cos2(6ij - A0) + 2) ) (24) 



Introducing (3j — \ cos{6j — A9)\, this can be written as 

2 / " 

(121 = -A9Hl{x) 6 + ^ (|k|/3| - 2|k|/32 + |k| + 2t/34 - 4</3| + 2i + t^/j, 

V j=i 

- 2i2/3| + - + bl3j - 5/3|t + 5t/3j - 4/3| + 2) ^ (25) 

Now, we integrate in k over the circle of radius K„iax- This translates to inte- 
grating |k| from to Kmax and multiplying by 27r.. This gives us 



+ ^ma. + 4i/3| - 8i/3| + 4< + 2i2/5^. - 4i2/3| 



+ 2^2 - 10/3| + 10/3j - 10/3j^t + lOtPj - 8/3f + 4^ | (26) 
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And this becomes 



12 + 4n + 2nr + Ant + K^^axn 



J2 {(Kma. + 'it)f3* - 10(t + 



3 = 1 



2{2t^ + At + K^ax + 'i)l3^ + 2(i' + 5t + 5)/?^) 



(27) 



Breaking up the sum yields: 

2 Krn n.cr. 



-A6'''Af (a;)(12 + 4n + 2nr + Ant + K„,axn 



{Kmax + 4i) - 10{t + 1)^/3: 



i=i 



- 2(2^2 + 4t + Kmax + 4) ^ /3| + 2(i2 + 5t + 5)J2 Pi) (28) 

Now, we notice that (3j — \ cos(/)| for some (j) and that | cos^| < 1. This finahy 
aUows us to remove the j dependence of the terms inside the sum, and we obtain, 
after substituting A6' back in: 



< 



■ ^ ,„ Af (a;) i2nt^ + 9nt + nK^ax + 8n + 3) (29) 



□ 



2.4 Improving the quadrature 

The result we describe here depends on the following idea: VT+^ has a branch 
cut on the region The rational function approximation can be expanded 

as a sum of first order poles, as per (|lip . Integrating an analytic function against 
this sum of poles (around a contour encircling [—«,«] yields a sum of the form 
J2n Wnf{zn), which approximates the integral of /(z) vl + z^ around a contour 
encircling [— z,i]. In particular, this is a second order quadrature. 

A natural line of inquiry is to ask is whether higher order quadratures can 
be used, simply by discarding the rational function approximation, and merely 
choosing a sum of poles according to some appropriate quadrature rule. We 
conjecture that this can be done. 
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